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The non-radiative scattering amplitudes for electron-positron annihilation into quark and lepton 
pairs in the TeV energy range are calculated in the double-logarithmic approximation. The expres- 
sions for the amplitudes are obtained using infrared evolution equations with different cut-offs for 
virtual photons and for W and Z bosons, and compared with previous results obtained with an 
universal cut-off. 
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^ , I. INTRODUCTION 

Next future linear e + e~ colliders will be operating in a energy domain which is much higher than the electroweak 
bosons masses, so that the full knowledge of the scattering amplitudes for e + e~ annihilation into quark and lepton 

■ pairs will be needed. The forward-backward asymmetry for e + e~ annihilation into leptons or hadrons produced at 
' energies much greater than the W and Z boson masses has been recently considered in Kef. |l| , where the electroweak 

radiative corrections were calculated to all orders in the double- logarithmic approximation (DLA) . It was shown that 
the effect of the electroweak DL radiative corrections on the value of the forward-backward asymmetry is quite sizable 
and grows rapidly with the energy. As usual, the asymmetry is denned as the difference between the forward and 

■ the backward scattering amplitudes over the sum of them. These amplitudes were calculated in Ref. |l| in DLA, by 
| introducing and solving the Infrared Evolution Equations (IREE) . This method is a very simple and the most efficient 

instrument for performing all-orders double- logarithmic calculations (see Ref. and Refs. therein). In particular, 
when it was applied in Ref. [2| to calculate the electroweak Sudakov (infrared-divergent) logarithms, it led easily to 
the proof of the exponentiation of the Sudakov logarithms. At that moment this was in contradiction to the non- 
exponentiation claimed in Ref. Q and obtained by other means. This contradiction provoked a large discussion about 
Oh, the exponentiation. The exponentiation was confirmed eventually by the two-loop calculations in Refs. Q-0 and by 
summing up the higher loop DL contributions in Refs. and These Sudakov logarithms provide the whole set 
of DL contributions to the 2 — > 2 amplitudes only when the process is considered in the hard kinematic region where 
all the Mandelstam variables s,t,u are of the same order. On the other hand, when the kinematics of the 2 — > 2 
processes is of the Regge type, besides the Sudakov logarithms, another kind of DL contributions arises, coming from 
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ladder Feynman graphs. Accounting for those (infrared stable) contributions it leads, instead of simple exponentials, 

m 



X 

to much more complicated expressions for the scattering amplitudes. This was first shown in Ref. |10j| . where in the 
framework of pure QED, the scattering amplitudes for the forward and backward e + e~ — » fi + fjT annihilation were 
calculated in the Regge kinematics. One example of high-energy electroweak processes in the Regge kinematics was 
considered in Ref. 0, where the backward scattering amplitude was calculated, for the annihilation of a lepton pair 
with same helicities into another pair of leptons. More general calculations of the forward and backward electroweak 
scattering amplitudes were done in Ref 0. 

Bowever, both calculations in Refs. and were done under the assumption that the transverse momenta ka_ of 
the virtual photons and virtual W, Z -bosons were much greater than the masses of the weak bosons. In other words, 
the same infrared cut-off M in the transverse momentum space, was used for all virtual electroweak bosons, i.e., 

h± > M > M w « M z . (1) 

Obviously, while M is the natural infrared cut-off for the logarithmic contributions involving W, Z bosons, the 
cut-off for the photons can be chosen independently, in accord with the experimental resolution in a given observed 
process. Indeed the assumption qTp - although simplifying the calculations a lot, is unnecessary and an approach that 
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involves different cut-offs for photons and W, Z weak bosons would be more interesting and suitable for phenomeno- 
logical applications. This technique involving different cut-offs for photons and for W, Z bosons was applied in Ref. 0, 
for calculating the double-logarithmic contributions of soft virtual electroweak bosons (the Sudakov electroweak log- 
arithms) but not for the scattering amplitudes in the regions of Regge kinematics. 

In the present paper we generalize the results of Refs. |2( and yj, and obtain new double-logarithmic expressions 
for the 2^2- electroweak amplitudes in the forward and backward kinematics. These expressions involve therefore 
different infrared cut-offs for virtual photons and virtual weak bosons. Throughout the paper we assume that the 
photon cut-off, fi, and the W, Z boson cut-off, M, satisfy the relations 

M>M WtZ , fJ->m f (2) 

where my is the largest mass of the quarks or leptons involved in the process. Notice that the values of M and /x could 
be widely different. Let us remind that in order to study a scattering amplitude A(s, t) in the Regge kinematics s>-f 
(where s and t are the standard Mandelstam variables), it is convenient to represent A(s,t) in the following form: 
A(s,t) = A {+ ^(s,t) +A ( .-^(s,t), with A^(s,t) = (l/2)[A(s,i) ± A(-s,i)\ called the positive (negative) signature 
amplitudes. We shall consider only amplitudes with the positive signatures. The IREE for the negative signature 
electroweak amplitudes can be obtained in a similar way, see e.g. Ref. pj for more details. 

The paper is organized as follows: in Sect. 2 we define the kinematics and express the scattering amplitude for the 
e + e _ annihilation in terms of invariant amplitudes. In Sect. 3, we construct the evolution equations for the invariant 
amplitudes for the case when in the center mass (cm) frame, the scattering angles are very small. First, we obtain the 
IREE equations in the integral form and then we transform them in the simpler, differential form. These differential 
equations are solved in Sect. 4 and explicit expressions for the invariant amplitudes involving the Mellin integrals 
are obtained. In Sect. 5, we consider the case of large scattering angles, or when the Mandelstam variables s, t and 
u are all large. Sect. 6 deals with the expansion of the invariant amplitudes into the perturbative series in order 
to extract the first-loop and the second-loop contributions. Then we compare these contributions to the analogous 
terms obtained when one universal cut-off is used and study their difference. The effect of high-order contributions 
in the two approaches is further studied in Sect. 7 where the asymptotic expressions of the amplitudes are compared. 
Finally, Sect. 8 contains our concluding remarks. 



II. INVARIANT AMPLITUDES FOR THE ANNIHILATION PROCESSES 



Let us consider a general process where the lepton I (pi) and its anti-particle U{p2) annihilate into a quark or a 
lepton q k (p[) and its anti-particle tfc'(p' 2 ) (see Fig. I): 



^(Pi)^(p2)^/(p'i)^b' 2 ) . (3) 



Pi Pi 




Pi P2 

FIG. 1: Scattering amplitude of the annihilation of Eq. 3. 



For this process, the most complicate case occurs when both the initial and the final particles (anti-particles) are 
left-handed (right-handed). The scattering amplitudes for other helicities can be obtained easily from the formulae 
derived for this case. As there is no technical difference when considering the annihilation into quarks or leptons, 
we present parallel results for the annihilation into a quark-antiquark or a lepton-antilepton pair. According to our 
assumption, the initial lepton belongs to the weak isodoublet (v : e). The final lepton belongs to another doublet, e.g. 
(z^t, fi), and the final quarks are also from a doublet, e.g. (it, d). The antilepton and the antiquark belong to the 
charge conjugate doublets. Obviously, the scattering amplitude A for the annihilation can be written as follows: 



A = q k '{p\)qAp' 2 )At k l k {p 1 )l i {p 2 ) , 



(4) 
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where the SU(2) matrix amplitude A k , k has to be calculated. We will consider it in DLA. The DL contributions 
to A", k are different according to the kinematics of the process. The kinematics is defined by appropriate relations 
among the Mandclstam variables s, t, u, 



= (Pi+P 2 ) 2 , ^(pi-p'i) 2 , u=( Pl -p' 2 ) 2 . 



Throughout this paper we assume that yfs ^> Mw,z- 
The kinematical regime defined as 



(5) 
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is called the hard kinematics and corresponds to large cm scattering angles 9 = 9 Pl ~ 1 ■ Radiative corrections to 
the annihilation in this kinematics yield DL contributions. 

There are also two other Regge-type kinematical regimes where DL contributions appear. First, there is the 
configuration where 



-u > -t 



(7) 



We call it the t -kinematics. According to the terminology introduced in 10], it is the forward kinematics (with 
respect to the charge flow) for e + e~ — * and e + e~ — > dd. At the same time, it corresponds to the backward 

kinematics for e + e~ — > uu. In this kinematics, 9 <C 1. 

Second, there is the opposite kinematics where 9 ~ tt and therefore 



We define the configuration © as the u-kinematics. It corresponds to the forward scattering for 



(8) 
uu and 



the backward scattering for e + e 



MM 



dd. 



To simplify the calculations, it is convenient to introduce the projection operators (-Fj)fc'ft (i = 1>2,3,4), so that 
A k , k can be written in the following form: 



A" — 

A k'k — 



u{-p' 2 )lp,u{p' 1 )u(-p 2 )l> 1 u(p l ) 



(P 



j+Dkk'Aj+i 



(9) 



where j ' — 1 for the i-kinematics and J = 3 for the it-kinematics. The representation (|5J) reduces the calculation of the 
matrix amplitude A k \ k to the calculation of the invariant amplitudes Aj. The explicit expressions for the operators 
(P c )k>k ( c = 1) •■'4) can be taken from Ref. Q: 



(10) 
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According to the results of Ref. Q, the forward (Ap) and backward (As) amplitudes of the e + e annihilation into 
quarks are expressed through invariant amplitudes Aj as follows: 



A F (e + e 
A F (e + e~ 



uu) = {A 3 + At)/2, A B (e + e -> uu) = M, 
dd) = (Ai + A 2 )/2, A B (e + e- -> dd) = A 2 . 



and the annihilation into lcptons is expressed through the leptonic invariant amplitudes very similarly: 



A F {e+e- -> n+ f i-) = (A 1 + A 2 )/2, A B (e + e" 
A F (e + e- -» v^) = (A 3 + At)/2, A B (e + e 



! = A 4 



(11) 



(12) 



We have used the general notation Aj for the invariant amplitudes in Eqs. (|11II12|) and we will keep using this notation 
until Sect. 4. 

In order to calculate the amplitudes Aj to all orders in the electroweak couplings in the DLA, we construct and 
solve some infrared evolution equations (IREE). These equations describe the evolution of Aj, (j = 1,2,3,4) with 
respect to an infrared cut-off. We introduce two such cut-offs, /z and M. We presume that M « M z « M w and use 
this cut-off to regulate the DL contributions involving soft (almost on-shell) virtual W, Z -bosons. In order to regulate 
the IR divergences arising from soft photons we use the cut-off fi and we assume that /i w m q <C M where m q is the 



4 



maximal quark mass involved. Both cut-offs are introduced in the transverse momentum space (with respect to the 
plane formed by momenta of the initial lcptons) so that the transverse momenta ki of virtual photons obey 

kij_> /i , (13) 

while the momenta ki of virtual W, Z -bosons obey 

k l± > M . (14) 

Let us first consider Aj in the collinear kinematics where, in the cm frame, the produced quarks or leptons move 
very close to the e + e _ -beams. In order to fix such kinematics, we implement Eq. J7J) by the further restriction on t: 

s u > M 2 > ju 2 > -t (15) 

and similarly for Eq. JSJ) by 

s 1 > M 2 > n 2 > -u . (16) 

Basically in DLA, the invariant amplitudes Aj depend on s,u and t through logarithms. Under the restriction 
imposed by Eqs. (|15l then all Aj depend only on logarithms of s,M 2 ,fj, 2 in the collinear kinematics. It is 
convenient to represent Aj in the following form: 

A 3 {s^ 2 ,M 2 ) = £D W 2 ) + ^(s,M 2 ,M 2 ) , (17) 

where Ay ED ' (s, /i 2 ) accounts for QED DL contributions only, i.e. the contributions of Feynman graphs without 
virtual W, Z bosons. To calculate A^ ED \s, ^i 2 ) we use the cut-off /i, therefore the amplitudes A < ^ ED ' > do not 
depend on M. In contrast, the amplitudes A' j(s, /i 2 , M 2 ) depend on both cut-offs. These amplitudes account for DL 
contributions of the Feynman graphs, with one or more W, Z propagators. By technical reasons, it is convenient to 
introduce two auxiliary amplitudes. The first one, Af ED \s,M 2 ), is the same QED amplitude but with a cut-off M. 

The second auxiliary amplitude, Aj(s,M 2 ) accounts for all electroweak DL contributions and the cut-off M is used 
to regulate both the virtual photons and the weak bosons infrared divergences. 

Beyond the Born approximation, the invariant amplitudes we have introduced depend on logarithms, the arguments 
of which can be chosen as in the following parameterization: 

A (QED) = A (QED ){s ^ 2) = A (QED) (s/ ^ j£QED) = ^D)^^ = Af ED )( S /M 2 ), (18) 

Aj = Aj(s,M 2 ) = A 3 (s/M 2 ), A' 3 = A' 3 {s,[i 2 ,M 2 ) = A'jis/M 2 ,^) , 

with 

r] = ln(M 2 //i 2 ) . (19) 

Our aim is to calculate the amplitudes A'j, whereas the amplitudes A^ ED \ A^ ED ^ and Aj(s/M 2 ) are supposed to 

be known. The amplitudes Aj were introduced and calculated in Ref. 0. In order to define amplitudes A'j, Aj(s, M 2 ), 
the projection operators of Eq. (fTUf> have been used. The use of these operators is based on the fact that the 
SU(2) x £7(1) symmetry for the electroweak scattering amplitudes takes place at energies much higher than the weak 

mass scale M. On the contrary, the QED amplitudes A^ ED ^ and A^ ED ^ are not SU(2) invariant at any energy. 
Nevertheless, it is convenient to introduce "the QED invariant amplitudes" A^ ED ^ , A^ ED ^ by explicit calculation of 
the forward and backward QED scattering amplitudes. Then inverting Eq. we construct the amplitudes A^ ED ^ 
for e + e~- annihilation into quarks: 

A[ QED] = 2Af ED \e + e- - dd) - A^ ED] (e + e~ - uu), Af ED) = A^ ED) (e + e- uu), (20) 
4 QED) = 2Af ED \e+e- - uu) - Af ED \e + e~ - dd), Af ED) = A^ ED \e+e- -> dd) 

and inverting Eq. I|12H allows us to obtain A^ ED ^ for e + e~- annihilation into leptons: 

A[ QED) = 2Af ED \e+e- -> Af ED) = 0, Af ED) = -A[ QED) = A { § ED \e + e- ^ + ,r). (21) 
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III. EVOLUTION EQUATIONS FOR AMPLITUDES A j IN THE COLLINEAR KINEMATICS 



We would like to discuss now the IREE for the amplitudes introduced earlier. The basic idea for constructing 
infrared evolution equations for the scattering amplitudes consists in introducing the infrared cut-offs in the transverse 
momentum space and evolving the scattering amplitudes with respect to them. This method does not involve analyzing 
the DL contributions of specific Feynman graphs but is based on quite general conceptions such as the analyticity of 
the scattering amplitudes and the dispersion relations which guarantees its applicability to a wide class of problems 
(see e.g. Ref. |2( and Refs. therein). The essence of the method is the factorizing the DL contributions of virtual 
particles with the minimal transverse momenta. The IREE with two cut-offs for the electroweak amplitudes in the 
hard kinematics 10 were obtained in Ref. . In the present section we construct the IREE for the 2 — ► 2 - electroweak 
amplitudes in the Regge kinematics. According to Eqs. (|13l I14|) . we use two different cut-offs for the virtual photons 
and for the weak bosons. The amplitude Aj is in the lhs of such an equation. The rhs contains several terms. In the 
first place, there is the Born amplitude Bj. In order to obtain the other terms in the rhs, we use the fact that the DL 
contributions of virtual particles with minimal transverse momenta (= k±) can be factorized. Furthermore, this k± 
acts as a new cut-off for the other virtual momenta. The virtual particle with k± (we call such a particle the softest 
one) can be either an electroweak bosons or a fermion. Let us suppose first that the softest particle is an electroweak 
boson. 

In this case, in the Feynman gauge, DL contributions come from the graphs where the softest propagator is attached 
to the external lines in every possible way whereas k± acts as a new cut-off for the blobs as shown in Fig. 2. When 



W,Z 




Symmetrical 
graphs 



FIG. 2: Softest boson contributions to IREE to A-j 



the softest electroweak boson is a photon, the integration region over k± is /i 2 <C k'j_ <C s and its contribution, Gj, to 
the rhs of the IREE is: 



=-^5 7) ( f^Hs/kl)Af ED \s,kl)+J^ ^Hs/kDA'^k^NP 

f ^M»/*i)4(«.*L*i)) , 

Jm* k± ' 



(22) 



where 



1(7) 2 • 2 a (Y2-Y1) 2 ( 7 ) 



7\V- 



g sin ow 



1 , 0-2 -vi 



(23) 



1,(7) _ „2 • 2 a ( F 2 + Y1) 2 ,(7) _ . 2 ,;,.2 , 



1W 



V? 1 = g z sin' 



1 , (Y 2 +Y 1 



We have used the standard notations in Eq. l|2"5|) : g, g' are the Standard Model couplings, Y\ (Y 2 ) is the hypercharge 
of the initial (final) fermions and 6\y is the Weinberg angle. The logarithmic factors in the integrands of Eq. 122|) 
correspond to the integration in the longitudinal momentum space. The amplitude A' in the last integral of Eq. 122|) 

does not depend on [i because k\ > M 2 . Therefore it can be expressed in terms of Aj(s, k 2 ^) and Ay ED '(s, k\) : 

A' j (s,kl,kl)=A j (s,ki)-Af ED \s,kl) . (24) 

When the softest boson is either a Z or a W, its DL contribution can be factorized in the region M 2 <§; k\ <C s. 
This yields: 

Gf Z) - ~^r ] SI If Hs/kDMs/kl) , (25) 



6 



with 



and the factors bj can be taken from Ref. [lj: 



= h - 6J (26) 



g (Yi — Y 2 ) 2 Sg 2 + g' 2 (Y 1 -Y 2 ) 2 _ 
61 = 4 &2 = ^ , (27) 

</ 2 (y 1 + y 2 ) 2 8, 9 2 + .9' 2 (y 1 + y 2 ) 2 

°3 = : , 04 = : ■ 



In Eq. I|25|) we have used the fact that the W and the Z bosons cannot be the softest particles for the amplitudes 
j^(QED) g j nce ^ ne integrations over the softest transverse momenta in Aj® ED ^ can go down to jj,, by definition. The 
sum of Eqs. I|22[) and (|25|l , Gj can be written in the more convenient way: 

G,(,s, M 2 , M 2 ) = Gf ( S ,/i 2 ,M 2 ) +Gf Z Hs,M 2 ) = Gf ED \s^ 2 ) - Gf ED \s,M 2 ) + Gj(s, M 2 ) + G'j(s, fi 2 ,M 2 ) , 

(28) 

where 

Gf ED) = -±b™ f^Hs/kl)Af ED \s/kl) , G f ^ = _ JjftW £ ^ln (s/fc i ) lf ^)( S /fci) , (29) 
6j = ~8^' H ln(.s/fcl)i,(.s/fci) , G'j = - ^ ln( S /fc 2 )^( S /fc 2 , M 2 /fc 2 ) . 

Eqs. I|28l I29|l account for DL contributions when the softest particle is an electroweak boson. However, the softest 
particle can also be a virtual fermion. In this case, DL contributions from the integration over the momentum fc of 
the softest fermion arise from the diagram shown in Fig. 3 where the amplitudes Aj are factorized into two on-shell 
amplitudes in the i-channel. We denote this contribution by Qj(s,fi 2 ,M 2 ). 

h 





Pi k -pi 

FIG. 3: Softest fermion contribution. 



The analytic expression for Qj is rather cumbersome. However it looks simpler when the Sudakov parameterization 
is introduced for the softest quark momentum k (with p\ and p 2 being the initial lepton momenta). 

k = ap 2 + Pl + k ± . (30) 

After simplifying the spin structure, we obtain 

Q j (s, fJ , 2 ,M 2 ) = c j j[dkl J ^ J {s J k \ ) 2 A i ^ k l ■ M *) A i k l ■ M ') - ( 31 ) 

where 

ci = c 2 = -c 3 = -c 4 = . (32) 
Similarly to Eq. Ij28(l . Qj of Eq. H31(l can be divided into the following simple contributions: 

Q 3 = Qf ED) -Qf ED) +Q ]+ Q'^ (33) 
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where 



er D V» 2 > - =, /; *i / 1 f p^4 8 * D Wi>4 8 * D > , (34) 
^"■c/iO = * jT *i / f f 5^5,^ • (») 

4(,,/A/ 2 ) = c, !' dk\ f B ,- IjtWfci ) , (36) 



M 2 J a (3 {sa(3 — k 2 ) 2 



and 



Q' j (s/M\r,) = c j j 2 dkl J ^ ^J^^ ^A^isa/kDA'jisP/k^Mykl) (37) 

+^' J (sa/fci,M 2 /fci)A' 3 (s/3/fci,Af 2 /fci)) . 
Now we are able to write the IREE for amplitudes Aj . The general form it given by: 

Aj = Bj + Gj + Qj . (38) 
Then using Eqs. I|28|) and (|33|l we can rewrite it as 

A'j + Af ED) = Bf ED) Bf ED) + Bj + Gf ED) - Gf ED) + Gj + G'j + Qf ED) Qf ED) + Qj + Q'j . (39) 
Let us notice that A^ ED \s/ fi 2 ) obeys the equation 

A (QED) = B (QED) + G (QED) + Q (QED) (4Q) 

and therefore A^ E cancels out in Eq. (|39")l . Also, the auxiliary amplitudes Aj and Aj ED \ obey similar equations: 

A (QED) = g(QED) + q(QED) + &{ QED) ^ = £ . + £ . + £ . (41) 

The solutions to Eqs. I|4(JI I41fl arc known. With the notations that we have used they can be taken directly from 
Ref. 0. Hence, we are left with the only unknown amplitude A'j in Eq. (|39|l . Using Eqs. (|40l I41f> . we arrive at the 
IREE for A'j, namely: 

A'j{s/M\ r,) = Aj(s/M 2 ) - Af ED \s/M 2 ) + G'j{s/M\r,) + Q'j(s/M 2 ,r,) . (42) 

In order to solve Eq. ()42[| . it is more convenient to use the Sommerfeld- Watson transform. As long as one considers 
the positive signature amplitudes, this transform formally coincides with the Mellin transform. It is convenient to use 
different forms of this transform for the invariant amplitudes we consider: 



7 CO 

too 



,( -' M3 >= /"s(^)"«">- (45) 



— too 

100 cL: 



A'j{s/M\ v ) = I ^(w) Fj(u,r)). (46) 
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Combining Eqs. I|43(l to (|46|l with Eq. I|42[) we arrive at the following equation for the Mellin amplitude Fj(u, if): 




where rf = ln(M 2 /fcjJ. Differentiating Eq. I|47|l with respect to /i 2 leads to the homogeneous partial differential 
equation for the on-shell amplitude Fj(ui,r)): 

^-5^(-| + "') + ^» + f)' < 48 » 

where we have used the fact that ln(s//j 2 ), in Eq. I|47|l. can be rewritten as ln(s/M 2 )+r] and that ln(s/M 2 ) corresponds 
to —d/duj . 



IV. SOLUTIONS TO THE EVOLUTION EQUATIONS FOR COLLINEAR KINEMATICS 

Let us consider first the particular case when b^ = 0. ft contributes to the forward leptonic, e + e~ — > fi + 
annihilation and corresponds, in our notations, to the option 



Y 1 =Y 2 = -1. (49) 

the forward e + e~ — > dd annihilation, though here Y\ = —1, Y<i = 

1/3 and therefore b^ 1 ' 7^ 0. In order to avoid confusion between these cases, we change our notations, denoting 

$1 = Fi, <f>i = fi and <^ 0) = / x (0) when Y 1 = Y 2 = 
discuss the annihilation into leptons. Then we deno 
the particular case (|49|l obeys the Riccati equation 

drj 



Let us notice that Aj with J = 1 contributes also to the forward e + e — > dd annihilation, though here Y\ = —1,12 
h h) _ 

$1 = Fi, 4>i = fi and </> 1 °' ) = f[ ^ when Y\ = Y2 = — 1. We will also use notations $2,3,4 instead of ^2,3,4 when we 
discuss the annihilation into leptons. Then we denote c = Ci = l/(87r 2 ). Therefore, the lepton amplitude $1(^,77) for 

quatii 

^=cM<V)* 1 + *0 , (50) 



with the general solution 



$1 C^-eNi'W (51) 
In order to specify C, we use the matching (see Eq. (14711 ) 

$i = 0i(w)-^i o) (u;), (52) 

when 7/ = 0, arriving immediately at 

1 ^)+^-(^-^))^% lMJ 
and therefore to the following expression for the invariant amplitude L\= A\ when Y± = Y 2 = — 1: 

du (s y ao), . r ^ / s f , 

Li - 2^ l^J ^ H + 2^ feJ #r^^v^ ■ ( } 

Obviously, when [i — > M, Eqs. I|54|l converges to the same amplitude obtained with using only one cut-off. Indeed, 
substituting \i — M and r\ = leads to 
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According to Eqs. (|21l I43f) . the QED amplitude is easily expressed in terms of Mellin amplitude cj^p 1 for the 
forward e + e~ — * annihilation: 

4 0) = 2^ 0) ■ (56) 
The expression for (jr F can be taken from Refs. [I3,|2 and 0: 



4°)=4 7 r 2 (o;- v /^- x 2) ) (57) 

with 

Xo - 2a/7r. (58) 
On the other hand, the amplitudes <f>j were calculated in Ref. p|. In particular, 



= Att 2 (lo~ v^^X 2 ) , (59) 
where % 2 is expressed through the electroweak couplings g and </': 

x 2 = [3g 2 + g' 2 }/(8n 2 ) . (60) 



Next, let us solve Eq. 1(48(1 for the general case of non-zero factor bj ■ Then, this equation describes the backward 
e + e~ annihilation into a lepton pair (e.g. /i + /i~) and also the forward and backward annihilation into quarks. Eq. 148(1 
looks simpler when uj, rj are replaced by new variables 

x = w/\j, y = \jT], (61) 



with Xj = yb^/(8ir 2 ). Changing to the new variables, we arrive again at the Riccati equation: 

**L = {a-T)F j -2 qj ff>F i - qi Ff , (62) 



where a = (x + y)/2, r = (x — y)/2 and = Cj/Xj. 
The general solution to Eqs. i(62(l is 



C{a) +qjQj(a,' 



Fi - ^ r:w,- , (63) 



where C(cr) should be specified, 



P 3 .(cr, r) = exp (<tt - t 2 /2 - 2qj f 



dUf\Q) (64) 



and 

QMr) = f + dCPjiaX) ■ (65) 

The QED amplitudes ff ] can be obtained from the known expressions for the backward, and forward /j ** 
QED scattering amplitudes: 

f£\x) = {Airae q /pf)d\u{e x2 ' A D {0) {x))/dx , (66) 

where D p are the Parabolic cylinder functions with pg = —2e q /(l + e q ) 2 and e q = 1 for the annihilation into muons, 
e q = 1/3 (2/3) for the annihilation into d (u)- quarks. Similarly, the QED forward scattering amplitudes for the 
annihilation into quarks are 

f$\x) = {AKae q lpf)d\n{e x2 ' A D< a) {x))/dx , (67) 
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with pg — 2e q /(l — e q ) 2 . Let us stress that the forward amplitudes for the annihilation into leptons are given by 

Eq. The amplitude fp B was obtained first in Ref. 0] for the backward scattering in QED. Obviously, the only 

difference between the formulae for fj(x) and ff\x) is in the different factors dj, Pj and Xj. We can specify C(a), 
using the matching 

F>) = />)-/j 0) M , (68) 
when T} = 0. The invariant amplitudes fj were calculated in Ref. pj: 



/<M = 2 f"°(^W) ^W). (69) 



Using Eq. we are led to 



(f J (x + y)-f^(x + y))P(a,r) 



* - -, : .. : , ', ^ (70) 



P(a, a) - + y) - ff\x + y)){Q{a, a) - Q(a, r)) 



and finally to 



/ -.x - " ' ■ M 

/-x 7 . . f /'.I , J_ ,,i _ 



It is easy to check that when p — M, Aj(s/M 2 ,rj) coincides with the amplitude Aj(s,M 2 ) obtained in Ref. Q. 

Eqs. I|54l I71|) describe all invariant amplitudes for e + e~ -annihilation into a quark or a lepton pair in the collinear 
kinematics H15II16JI . 



V. SCATTERING AMPLITUDES AT LARGE VALUES OF t AND u 

In this section we calculate the scattering amplitudes A when the restriction of Eqs. (|15l I16|) for the kinematical 
configurations (Q |SJl are replaced by 

s > Af 2 > -t > n 2 (72) 

and 

s > M 2 >-u^>p 2 . (73) 

In this kinematical regions it is more convenient to study the scattering amplitudes A directly, rather than using 
the invariant amplitudes Aj . In order to unify the discussion for both kinematics (|72l 17311 , let us introduce 

k = — t , (74) 

when ij?^ is considered and 

« = — u (75) 

for the other case (J73J). Using this notation, the same parameterization A = A(s, p 2 , M 2 , ft) holds for both kinemat- 
ics <|72l I73|) . Let us discuss now the evolution equations for A. As in the previous case, it is convenient to consider 
separately the purely QED part, A^ ED ^ and the mixed part, A': 

A(s, ft, p 2 , M 2 ) = A^ QED '>(s 1 ft, p 2 ) + A'(s, ft, m 2 , M 2 ). (76) 

Generalizing Eq. i|18|) . we can parameterize them as follows: 

A^ ED \s,K,p 2 ) = A^ ED \s/p 2 , K /p 2 ), A'( s , K ,p 2 ,M 2 )=A'( s /M 2 ,s/p 2 ,K/[i 2 7 M 2 /p 2 ). (77) 
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In order to construct the IREE for A^ ED ^ and A', we should consider again all options for the softest virtual 
particles. The Born terms for the configurations l|72|l and (|73H do not depend on p 2 and vanish after differentiating 
on p. The same is true for the softest quark contributions. Indeed, the softest fermion pair yields DL contributions 
in the integration region k\ 3> k, which is unrelated to p. Hence, we are left with the only option for the softest 
particle to be an electroweak boson. The factorization region for this kinematics is 

p 2 < k\ < k . (78) 

Obviously, only virtual photons can be factorized in this factorization region, which leads to a simple IREE: 

^ + ^ = -^ V + ^ M «->, ™ + ^ + §£ = -^V + ^< (79) 

where we have denoted p = \n(s/p 2 ), z = \ii(k/p 2 ), t]' = 111(77) = hi(M 2 //z 2 ) and A = a/2-K. The factors and h^ 1 ' 
are: 

feW = eie ' 1+ e 2 e / 2, 6 (7) = e x e 2 + e\e' 2 - e 2 e\ - e x e>i (80) 

for the case {72J, and 

= — e 2 e'i + eie'2, b' 7 ' = eie 2 + e'ie' 2 + e 2 e' 2 + eie'i (81) 

for the other case (|75|) . 

The notations ei,e'i in Eqs. I|80l I81() stand for the absolute values of the electric charges. They correspond to the 
notations of the external particle momenta introduced in Fig. 1. The terms proportional to hfr' in Eq. (|79|) correspond 
to the Feynman graphs where the softest photons propagate in the K-channels. Let us notice that for any kinematics 
we consider it holds 

b^ + h™ = {l/2)[e\ + el + e'l + e'l] (82) 

due to the electric charge conservation. 

In order to solve Eq. I|79J) . we use the matching with the amplitude A(s, p 2 , M 2 ) for the same process, however in 
the collinear kinematics: 

A^ ED \ s ,p 2 7 K,M 2 )=A^ ED \s,p 2 ), A'(s, K ,p 2 7 M 2 )=A'(s,p 2 ,M 2 ) , (83) 

when k — p 2 . The solution to Eq. (|79H is 

A( QED) = lp (QED) ip _ z)e -X b ^^ / 2-Xh^^/^ A > = ^ {p - Zjf] > z)e ~^ P*/*-^ I* . (84 ) 

Using the matching of Eq. (|83[l allows to specify tp and ip^ ED \ After that we obtain: 

a (qed) = S 'A^ ed \s/k), A' = S'A'(s/M 2 ,M 2 /k) , (85) 

where 

g> = e -M™pz+\{b™-h™)*>/7 (86) 

We did not change s/M 2 to s/k in Eq. Ij85(l because M 2 ^> k. It is convenient to absorb the term — Xb^pz into the 
amplitudes A^ ED ^ and A'. Introducing, instead of uj, the new Mellin variable I = uj + Xb^z (see Ref. jj for details), 
we rewrite Eq. <|85|l as follows (for the sake of simplicity we keep the same notations for these new amplitudes A^ ED ^ 
and A'): 

a {qed) =S (A( QED \ s / K ) + A'{s,K, i j?,M 2 )) (87) 

with S being the Sudakov form factor for the case under discussion. S includes the softest, infrared divergent DL 
contributions. When the photon infrared cut-off p is assumed to be greater than the masses of the involved fermions, 
this form factor is: 

S = exp ( - ^(b^ + fcW) ln 2 ( K / M 2 )) . (88) 
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However, in the case of e + e annihilation into quarks (muons), if the cut-off fj, is chosen to be very small, less than 
the electron mass, m e the exponent in Eq. (|88|) should be changed to : 

S = exp ( - + h^)(ln 2 ( K /n 2 ) - ln 2 (m 2 J^ 2 ) - ln 2 (m 2 / M 2 )) , (89) 

where m is the mass of the produced quark or lepton (cf . Ref . |ll| ) . 

If m > fi > m ei the last term in the exponent of Eq. (|89|l is absent. The kinematics with larger values of /?, e.g. 
s>k» M 2 , can be studied similarly, although it is more convenient to use the invariant amplitudes Aj. The result 



is 



A^SMs/k) , (90) 



where 

Sj = exp 



((^ 7) + ^. T) )(ln 2 ( K //i 2 ) - ln 2 (m 2 / M 2 ) - ln 2 (m 2 /fi 2 ) + (b, - ft™ + hj - h™) ln 2 ( K /A/ 2 ))] (91) 



and A(s/M 2 ) is the scattering amplitude of the same process in the limit of collinear kinematics and using a single 
cut-off M. These amplitudes were defined in Sect. 2. The factors hj given below were calculated in Ref 

h x = g 2 (3 + tan 2 6 w Y x Y 2 )/2, h 2 = g 2 (-l + tan 2 6 w Y 1 Y 2 )/2 , (92) 
h 3 = g 2 (3 - tan 2 e w Y l Y 2 )/2, h 4 = g 2 (-l - tan 2 6 w Y 1 Y 2 )/2 . 

The form factors S, Sj include the soft DL contributions, with the cm energies of virtual particles ranging from \j? 
to k. Due to gauge invariance, the sums bj + and bj + hj do not depend on j and Sj is actually the same for 

every invariant amplitude contributing to A^, k in the forward (backward) kinematics (see Ref. [lj). Obviously, in the 
case of the hard kinematics where (see Eq. (0) s ~ — u ~ — t, i.e. s ~ K, ladder graphs do not yield DL contributions. 
The easiest way to see this, is to notice that the factor (s/k) u in the the Mellin integrals l|45|) for amplitudes Aj does 
not depend on s in the hard kinematics, therefore all Mellin integrals do not depend on s. So, the only source of DL 
terms in this kinematics is given by the Sudakov form factor Sj given by Eq. (|91|l . Therefore, we easily arrive at the 
known result 

A k'k = B k'k s j ■ (93) 

B % y k in Eq. H93[) stands for the Born terms. The electroweak Sudakov form factor l|91|l with two infrared cut-offs was 
obtained in Ref. @- 

VI. FORWARD e + e~ ANNIHILATION INTO LEPTONS 

Eqs. Q11II12[I53|) and J7TJ give the explicit expressions for the scattering amplitudes of e + e -annihilation into quarks 
and leptons in the collinear kinematics. These expressions resume the DL contributions to all orders in the electroweak 
couplings and operate with two infrared cut-offs. In order to estimate the impact of the two-cuts approach, we compare 
these results to the formulae for the same scattering amplitudes obtained in Ref. Q where one universal cut-off M 
was used. We focus on the particular case of the scattering amplitudes for the forward e + e~ annihilation into leptons 
and restrict ourselves, for the sake of simplicity, to the collinear kinematics of Eq. (|15fl . Other amplitudes, and other 
kinematics can be considered in a very similar way. Eqs. 1121 1541 and (|71|) show that the scattering amplitude 
of the forward e + e~ into [i~ fi + is 



- 



1 r°° ,j_y fc{x + y)P 2 {a,T) 

2j_ loo 2 m \M 2 J P 2 (a,o)- ( j )2 {x + y )[Q 2 (a, ( j)~Q 2 ( ( j,T)\ ' 

The first integral in this equation accounts for purely QED double-logarithmic contributions and depends on the 
QED cut-off /i whereas the next integrals sum up mixed QED and weak double-logarithmic terms and depend on 
both n and M. The first and the second integrals in Eq. (|94|l grow with s whilst the last integral rapidly falls 
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when s increases. The point is that this term actually is the amplitude for the backward annihilation into muon 
neutrinos. It is easy to check that the QED amplitudes <j)p^ vanish when /i = M and the total integrand contains only 
[4>i (id) + (j>2(u))]/2. In contrast to Eq. (|94|l . purely QED contributions are absent in formulae for e + e~ annihilation 
into neutrinos. For example, the scattering amplitude L^) of the forward e + e~ — * v^v^ -annihilation in the collinear 
kinematics is 



r w _ 1 r° d ^ ( S V 



<j> 3 (x + y)P 3 (a,r) | 



P a (a, a) -(j) 3 (x + y)[Qs((7,a) - Q 3 (a,r)} 
4>a(x + y)Pi(^, r) 



, ^ 27 r 2 (^ + ^ + ^ + ...) , (96) 



2tt 2 (- 



P 4 (<j, a) - <j) A (x + y)[Q4,(cr,a) - Q 4 (a,r)] 

Similarly to Eq. 1(94(1 . the integrand in Eq. I|95() is equal to [03 (id) + 04(w)]/2 when /i — M. Although formally 
Eqs. I|94l 195(1 correspond to the exclusive e + e~ annihilation into two leptons, actually these expressions also describe 
the inclusive processes when the emission of photons with cm energies < [i is accounted for. 

Let us study the impact of our two-cut-offs approach on the scattering amplitude L^f' of Eq. 1(94(1 . As the last 
integral in Eq. 1(94(1 rapidly falls with s, it is neglected in our estimates and we consider contributions of the first and 
the second integrals only. First we compare the one-loop and two- loop contributions. Such contributions can be easily 
obtained expanding the rhs of Eq. ((94(1 into a perturbative series. From Eqs. 1(57(1 and ((59(1 one obtains that 

\ u Aid 3 8 id 
id 4cj 3 8^ 5 "7 ' 

with xo> X defined in Eqs. ((581 lt)0|) . Substituting these series into the first and the second integrals of Eq. ((94|) and 
performing the integrations over u>, we arrive at 

L« = 7 { 1} ln 2 (s/ M 2 ) + 7 < 1} ln(s//i 2 ) ln(s/M 2 ) + 7 | X) ln 2 (s/M 2 ) (97) 

for the first-loop contribution to Lp and 

L (2) = 7 p ) ln 4 (s/ M 2 )+ 7 ( 2) ln 3 (s/ M 2 )ln(s/M 2 )+ (98) 
7l 2) ln 2 (s/ M 2 ) ln 2 (s/M 2 ) + 7 f } \n(s/^ 2 ) ln 3 (s/M 2 ) + 7 < 2) ln 4 (s/M 2 ) 

(k) 

for the second-loop contribution. The coefficients 7^- are given below: 

^!V tf L ^;'»') ,^- ^;'»'l (99, 

(2) _ 7T 2 XO 6 (2) _ (2) _ (X 4 -4 X Q 4 ) 

7i --gg-.Ta -0,73 - 32 > 

(2) = ^ 2 (x 6 -6x 2 xo 4 + 2xo 6 ) (2) = ^ (3x 6 -24 X 2 xo 4 + 14xo 6 ) 
7 4 24 1 7 5 lg2 

Let us compare the above results with those obtained with one universal cut-off M only. We introduce the notation 
L(s/M 2 ) for amplitude when one cut-off M is used. The ratio R 1 = L 1 (s,fx,M)/L^(s,M) of the first loop 
contributions to the amplitudes and L is 

r(l) 

i?« = £ (100) 

^\n 2 (s/M 2 ) V ' 

where 7 1 = 7r 2 x 4 /8. Similarly the ratio R^ of the second-loop contributions is 

r(2) 

i? (2) = r - (101) 

^\n A (s/M 2 ) V 7 

where 7 2 = 7r 2 x 6 /64. Eqs. 1(1001 HOlf show explicitly that the difference between the one cut-off amplitude L and the 
two cut-off amplitude Lp^ grows with the order of the perturbative expansion, though rapidly decreasing with s. 

We can expect therefore that a sizable difference between and L when all orders of the perturbative series are 
resumed. 
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12 3 4 

FIG. 4: Dependence of i?' 1 ' on s for different values of /i(GeV). 




12 3 4 

FIG. 5: Dependence of R> 2 ' on s for different values of /x(GeV). 



VII. ASYMPTOTICS OF THE FORWARD SCATTERING AMPLITUDE FOR e+e^ ANNIHILATION 

INTO fi~ . 

In order to estimate the effect of higher order DL contributions on the difference between the one-cut-off and two- 
cut-off amplitudes, it is convenient to compare their high-energy asymptotics. For the sake of simplicity, we present 
below such asymptotical estimates for the amplitude L F of the forward e + e~ annihilation into fi + in the collinear 
kinematics (|15H . Calculations for the other amplitudes l|71f) can be done in a similar way. As well-known, the leading 
contribution to the asymptotic behavior is L F ~ s"°, with uq being the rightmost singularity of the amplitude L F . 

This amplitude contains the amplitudes <j>f\ and and therefore also their singularities. Eqs. (|59l I57|) show that 

the singularities of both 4>i and (jJf^ are the square root branching points. The rightmost singularity of (f>^ is xo an d 
the rightmost singularity of <j>i is \. They are defined in Eqs. 1581 16U|) . Obviously, 

0i O) (xo)=4^ 2 Xo, ^ 0) (x) = 4tt 2 ( X - tJx 2 - Xl) = 4vr 2 (x - x') , Mx) = ^X ■ (102) 
Combining Eqs. I|94|) and I|1U2|) and neglecting the last integral in Eq. I|94|l . we obtain the asymptotic formula for 
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the forward leptonic invariant amplitude A: 



Xo + 4tt' 



/ s y 2( X -x / )(2x , -x)e 2T)(x ' x ' ) 
\M 2 J 3x-2 X ' - (2x' - x)e 2r ' (x " x,) ' 



(103) 



The first term in Eq. (|l(J3f) represents the asymptotic contribution of the QED Feynman graphs, the second term 
the mixing of QED and weak DL contributions. On the other hand, when the one-cut-off approach is used, the new 
amplitude L F asymptotically behaves as: 



Then defining Z(s,rj), as: 



it is easy to see that 



Z(S,T)) 



F 



fj_y x+xo ^xo e r, 



\M 2 ) 



4(X-X / )(2X , -X)e 2r ' (x - X ' ) 
X[3x-2x'~ (2 X ' - x)e 2 "^')] 



(104) 



(105) 



(106) 



As Xo < Xi Z( s ) f &us when s grows. So, the one-cut-off and the two-cut-off approach lead to the same asymptotics, 
although at very high energies, say y/s > 10 6 TeV. At lower energies, accounting for Z, the amplitude is increased 
by a factor of order 2. On the other hand, Z strongly depends on the ratio M/n, which, of course, is related to the 
actual phenomenological conditions. To illustrate this dependence, we take M = 100 GeV and choose different values 
for /X, ranging from 0.1 to 1 GeV. Then in Fig. 6 we plot Z(s,fi) for fj, = 1 GeV and /i = 0.5 GeV. This shows that 
the variation is approximately 1.5 at energies in the interval from 0.5 to 5 TeV. 



1.5 

1.4; 

1.3 
1.2 
1.1 

1 ■ 
0.9 



fj, = 0.5 



M = 1 



2 

V^(TeV) 



FIG. 6: Dependence of Z on s for different values of \x (GeV). 



It is also interesting to estimate the difference between the purely QED asymptotics of L F (the first term in the 
rhs of Eq. I|103fl ') and the full electroweak asymptotics. To this aim, we introduce Aew'- 

L^ = (L^) {QED \l + A EW ). (107) 
From Eq. (|103fl we immediately get the following asymptotic behavior for Aew : 

a ( s Y- xo 2(x-x , )(2x / "X)e 2 " (x - x) nns , 

Aew ~ lil^J 3x-2x'-(2x'-x)e 2 ^') ( } 

As x > Xo; Aew grows with s, as shown in Fig. 7, Therefore the weak interactions contribution is approximately of 
the same size of the QED contribution, and their ratio rapidly increases as /x decreases. 
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12 3 4 

FIG. 7: Dependence of Aew on s for different values of [i (GeV). 



VIII. SUMMARY AND OUTLOOK 

Next future linear e + e~ colliders will be operating in a energy domain which is much higher than the electroweak 
bosons masses, so that the full knowledge of the scattering amplitudes for e + e~ annihilation into fermion pairs will 
be needed. In the present paper we have considered the high-energy non-radiative scattering amplitudes for e + e~ 
annihilation into leptons and quarks in the Regge kinematics Q and ijHJ). We have calculated these amplitudes in the 
DLA, using a cut-off M, with M > Mz ~ M^, for the transverse momenta of virtual weak bosons and an infrared cut- 
off /i for regulating DL contributions of virtual soft photons. We have obtained explicit expressions l|53l I71fl for these 
amplitudes in the collinear kinematics I|15lll6fl and Eqs. (|87II90|I for the configuration where all Mandelstam variables 
are large. The basic structure of the expressions in the limit of collinear kinematics is quite clear. They consist of two 
terms: the first term presents the purely QED contribution, i.e. the one with virtual photon exchanges only, whereas 
the next term describe the combined effect of all electroweak boson exchanges. Obviously, in the limit when the cut-off 
fi — > M, our expressions for the scattering amplitude converge to the much simpler expressions obtained in Ref Q 
with one universal cut-off for all electroweak bosons. In order to calculate the electroweak scattering amplitudes, we 
derived and solved infrared equations for the evolution of the amplitudes with respect to the cut-offs M and fj,. 

In order to illustrate the difference between the two methods, we have considered in more detail the scattering 
amplitude Lp of the forward e + e _ annihilation into and studied the ratios of the results obtained in the 

two approaches, first in one- and two-loop approximation and then to all orders to DLA. The ratios of the first- 
and second-loop DL results are plotted in Figs. 4 and 5. The total effect of higher-loop contributions is estimated 
comparing the asymptotic behaviors of the amplitudes. This is shown in Fig. 6. The effect of all electroweak DL 
corrections compared the QED ones is plotted in Fig. 7. It follows that accounting for all electroweak radiative 
corrections Ljf* increases by up to factor of 2.5 at y/s < 1 TeV, depending on the value of Mj ^i. In formulae for the 
2^2- electroweak cross sections, one can put M = Mw ~ Mz whereas the value of fj, is quite arbitrary. However 
it vanishes, when these expressions are combined with cross sections of the radiative 2 — > 2 + X processes. 

In the present paper we have considered the most complicated case of both the initial electron and the final quark or 
lepton being heft-handed (and their antiparticles right-handed). Studying other combinations of the helicities of the 
initial and final particles can be done quite similarly. We intend to use the results obtained in the present paper for 
further studying the forward-backward asymmetry at TeV energies, by including also the real radiative contributions. 
Basically, the QCD radiative corrections can give a big impact on the amplitudes of e + e~ - annihilation into hadrons. 
However, the perturbative QCD corrections cancel out of the expressions for the forward-backward asymmetry (see 
Ref. 0) whereas the non-perturbative corrections describing hadronization of the produced qq - pairs can be accounted 
for in the same way as it was done in Ref. |l(. 
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